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Abstract. We report on deep BVR-imaging of the field of the nearby millisecond pulsar PSR J0030+0451 obtained 
with the ESO/VLT/FORS2. We do not detect any optical counterpart down to B > 27.3, V > 27.0 and R > 27.0 
in the immediate vicinity of the radio pulsar position. The closest detected sources are offset by > 3", and they 
are excluded as counterpart candidates by our astrometry. Using our upper limits in the optical, and including 
recent XMM-Newton X-ray data we show that any nonthermal power-law spectral component of neutron star 
magnetospheric origin, as suggested by the interpretation of X-ray data, must be suppressed by at least a factor 
of ~ 500 in the optical range. This either rules out the nonthermal interpretation or suggests a dramatic spectral 
break in the 0.003 — 0.1 keV range of the power- law spectrum. Such a situation has never been observed in 
the optical/X-ray spectral region of ordinary pulsars, and the origin of such a break is unclear. An alternative 
interpretation with a purely thermal X-ray spectrum is consistent with our optical upper limits. In this case the 
X-ray emission is dominated by hot polar caps of the pulsar. 
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1. Introduction 

Millisecond pulsars (hereafter MSPs) differ from ordi- 
nary radio pulsars by much shorter spin periods P, 
smaller period derivatives P, higher dynamical ages r, 
weaker magnetic fields B, and evol ution histories (see, 
e.g., recent review by Lorimer 2001 ). Contrary to ordi- 
nary pulsars, only 9 of 56 MSPs currently known in the 
Galactic disk and 25 o f 52 MSPs fou n d in globular cluster s 
are isolated objects ( Lorimer 200lt [Lorimer et al. 20021) ; 
Possenti et al. 2001 ). It is believed that the fast rotation 
of these neutron stars (NSs) was gained in the past by 
angular momentum transfer during mass accretion from 
a companion star ( Bhattacharya fc van den Heuvel 1991 ) . 



dozen radio MSPs have been detected in X-rays. It is re- 
markable that their efficiency in converting spindown en- 
ergy to X-ray luminosity is roughly the same as for or- 
dinary pulsars, Lx /E ~ 10^"^ ( Becker fc Triimper 1997 ; 
Becker et al. 2000| ) . This suggests that the emission mech- 
anisms responsible for the multi-wavelength radiation 
of MSPs and ordinary pulsars can be similar, and one 
could therefore expect to detect MSPs in other spec- 
tral ranges as well, as has been done for several ordi- 
nary pulsars. Detection of the first MSP in gamma-rays 



This was supported by the discoveries of three accretion- 
powered X-ray MSPs in low-mass X-ray binaries (e .g., 
SAX J1808.4-3658, see |Wijnands fc van der Klis"T99^ ). 

Despite these differences, the distribution of integrated 
radio luminosities, as well as the luminosity dependence 
on P, P, B, and spindown energy losses E, are apparently 
similar for these much older and low-magnetized NSs, and 



for ordinary pulsars (Kuzmin & Losovsky 2001). About a 
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(Kuiper et al. 2000) supports this idea. 

To our knowledge, there are still no reports on opti- 
cal detection of isolated MSPs. It is hardly possible to 
detect thermal emission from the entire surface of these 
old, 10^ — 10^" yr, and cold NSs. However, the spindown 
energy, expected to power the nonthermal emission of pul- 
sars, can be much higher for MSPs than for old ordi- 
nary pulsars, and may even rival that of young pulsars. 
Assuming the same efficiency of conversion of spindown 
energy to nonthermal optical luminosity as for ordinary 
pulsars, one can estimate that nearby MSPs may well be 
detectable in the optical with large telescopes. A prob- 
lem is, however, that most of the nearby and energetic 
MSPs are components of close binary systems where the 
companion is predominantly either a white dwarf or main 



sequence star (Lorimer 2001) which outshines the pulsar 



in the optical. Fortunately, there are at least nine solitary 
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Table 1. Parameters of PSR J0030+0451 (Lommen et al. 2000 and Becker et al. 200C, unless specified otherwise). 







Observed 








Derived 




p 

ms 


P 

^0-20 


«2000, (^jooo ^J.a, Mi 

mas yr^^ 


I, h"^ 


DM" 
cm~'^ pc 


r 
Gyr 


B E 
G erg 


pc 


4.865 


f.0±0.2 


00''30'"2r:432(9)'= 4.1 ±3.4 
04°51'39'.'67(2) -24.6 ±8.0 


113? 1 

-57?6 


4.3328(2) 


7.7 


2.2 X 10* 3.4 X 10^^ 


230 



" Coordinates are at the epoch of the VLT observations, MJD 52134 (Aug 13, 2001) Galactic coordinates 

' Updated values of the proper motion (A. Lommen 2001, private communications) Dispersion measure 



Numbers in parentheses are uncertainties referri ng to the last significan t digit quoted 



^ _DM-based distance (the new distance model by Cordes & Lazio (2002) places the pulsar to 317 pc) 



MSPs in Galactic disk whose companions are believed to 
have been either evaporated or ablated (see, e.g., Lommen 
et al. pOOOl). 



Table 2. Log of VLT observations of PSR J0030+0451 in 
the BVRs bands. Exposure time is 720 s for all frames. 



Here we report on deep BVR imaging of the field of one 
of these solitary millisecond pulsars, PSR J0030+0451. 
This pulsar was only recently discovered with the 
Arecibo telescope (Lommen et al. 2000] ), and soon there- 
after detected in X-rays during the final observations 
with the ROSAT/PSPC ( [Becker et al. 2000) ). Recently 
it was re-observed in X-rays with the XMM-Newton 



(Becker & Aschenbach 2002). This relatively nearby NS 
(see Table |l| for its parameters) is characterized by high 
X-ray fiux, about (4.3 — 6.8) x 10"^'^ erg cm~^ in 
the 0.1 — 2.4 keV band, and low interstellar absorption, 
Nn ^ 3 X 10^" cm~^, corresponding to a color excess 
E{B — V) ^ 0.06 mag. This makes it a promising can- 
didate for optical detection. In Sect. 2 we present the ob- 
servations and the data reduction. In Sect. 3 we discuss 
our results in the optical in conjunction with the available 
X-ray data. 



2. Observations and data reduction 

The field of PSR J0030±0451 was observed in service 
mode on July 26 and August 13, 2001, with the FOcal 
Reducer/low dispersion Spectrograph (F0RS2) on the 
ES0/VLT/UT2 telescope[|, with a pixel scale of 0'.'2. We 
used Bessel filters for B and V, and an ESO special fil- 
ter for R (Rspociai, henceforth called Rs)^- Unfortunately, 
some of the V images were corrupted by bad CCD columns 
near the expected position of the pulsar and were not used 
in the analysis. The images we used are listed in Table 

Bias subtraction and fiatfielding were performed in a 
standard way, and the reduced individual images were 
aligned using a set of bright, non-saturated field stars. 
Standard utilities from the NOAO I RAF package were then 
used to combine the images applying the averaged sigma 
clipping algorithm avsigclip with the scale parameter equal 
to none. The pulsar vicinity is shown for B, V and Rs 
bands in Fig. |l|. 

^ See http://www.eso.org/instruments/fors/ for details on 
the instrument. 
^ See http://www.eso.org/instruments/forsl/Filters/. 
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Band 


Time 


Airmass 


Seeing 




UT 




UT 




arcsec 


1 


26.07.01 


B 


09:40 


1.172 


0.6 


2 






09:53 


1.188 


0.5 


3 


13.08.01 


B 


04:55 


1.579 


1.1 


4 






05:09 


1.495 


1.1 


5 






06:34 


1.208 


0.8 


6 






07:56 


1.150 


0.7 


7 






08:09 


1.155 


0.9 


8 
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05:38 


1.358 


0.7 


9 






05:52 


1.311 


0.8 


10 






06:06 


1.269 


0.8 


11 






07:14 


1.159 


0.7 


12 






07:28 


1.151 


0.6 


13 






07:42 


1.149 


0.7 


14 






06:47 


1.187 


0.7 


15 






08:24 


1.167 


0.7 


16 






08:37 


1.181 


0.7 


17 






08:53 


1.205 


0.7 


18 






09:07 


1.230 


0.7 



The radio position of PSR J0030±0451 at the epoch of 
the VLT observations (for which 13.08.01 was adopted, see 
Table 0) was determined using recent radio ephemerides 
(A. Lommen 2001, private communication for additional 
Arecibo observations). Astrometrical referencing of our 
images was made with I RAF tasks ccmap/cctran using 
the positions of several dozens of reference stars from the 
USNO A-2.0 catalogue seen in the images and optimizing 
the astrometrical fit by removing step by step the refer- 
ence stars with the largest residuals. For the 5 most suit- 
able stars we finally got Icr rms-errors of O'.'OS and O'.'ll, 
and maximum residuals of 0'/20 and 0'.'24, in RA and Dec, 
respectively. Combining the rms-errors with the nominal 
USNO accuracy of 0'/24 and radio ephemeris uncertain- 
ties, we obtained the 3(j pulsar VLT/FORS position un- 
certainties 0'/79 and 0'.'88 in RA and Dec, respectively. The 
resulting 3ct error ellipse is marked in the smaller images 
in Fig. |l|. No reliable counterpart candidate was detected 
within, or close to the error ellipse of the expected pulsar 
position. We also double-checked the astrometry with 5 
stars from the GSC-II catalog. The obtained rms errors 
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Fig. 1. Left panel: 1' x 1' overview of the 
PSR J0030+0451 field extracted from the 6.'8 x 6.'8 
frame of the VLT/F0RS2 image in the Rs band. The 
central box is blown-up in the smaller panels. A nearby 
bright USNO star (object 1) and two other stellar 
objects (2 and 5 from Table ^ are marked. Smaller 
panels: 8" x 8" pulsar vicinity in the BVRg bands. 
3(7 uncertainty ellipses of the expected pulsar position 
at the epoch of the VLT observations are marked. No 
pulsar is detected in these frames. Two objects are 
marked "5" (referring to Table and "6" in the Rs 
image, and are discussed in the text. 



Table 3. Magnitudes of several stars in the PSR J0030+0451 field with coordinates and offsets from the pulsar position 
shown in the table. Stars 1, 2 and 5 are marked within the left frame in Fig. |l|, while 3 and 4 are outside this frame. 



No 


02000 


52000 


Offset 


B 


V 





1 


00:30:27.83 


+04:51:35.2 


7'.' 5 


19.20+0.01 






2 


00:30:29.25 


+04:51:18.3 


34'.' 6 


22.32+0.02 


20.92+0.01 


20.27+0.01 


3 


00:30:39.33 


+04:53:15.4 


3' 21'.' 9 


23.77+0.06 


22.08+0.01 


21.29+0.01 


4 


00:30:25.39 


+04:53:14.2 


1'39'.'3 


24.15+0.07 


22.40+0.02 


21.78+0.01 


5 


00:30:27.32 


+04:51:36.9 


3'.'2 


23.82+0.05 


23.69+0.06 


23.63+0.06 


PSR 


00:30:27.43 


+04:51:39.67 




>27.3 


>27.0 


>27.0 



were O'.'OS and C/.'ll in RA and Dec, respectively, and the 
expected pulsar position was moved O'.'OG west and 0'.'14 
north in respect to the USNO position. However, since the 
main source of errors is the catalog uncertainty, we accept 
the USNO results as more conservative estimate. 

For the photometric calibrations we used the pho- 
tometr ic standards from the PG1323-085 and SA109 
fields ( |Landolt 1992 ), observed at the second night of 
our observations, and average Paranal extinction coef- 
ficient^ We then derived 3tT detection limits as m = 
-2.5 log (SCTV^/Texp) + mo, where a is the standard de- 
viation of the flux in counts per pixel, T^^p is the exposure 
time, A is the area of an aperture (in pix^) with a radius 
of 1" (corresponding to ^ 83% of the flux in a PSF of 
our images), and mo is the photometric zeropoint, includ- 
ing corrections for atmospheric extinction. The limits are: 
B = 27.3, V = 27.0, and = 27.0. 

In Table | we list BVRs magnitudes for several of 
the objects in the PSR J0030+0451 field. Object 1 is the 
USNO star U0900_00118426 marked in the left panel of 
Fig. |l|. It is oversaturated in our V and Rs images. Object 
5 is the source seen at the bottom right of the blown-up 
images in Fig. |l]. This is the closest object to the radio 



position of the pulsar clearly detected in all our images. 
Although it is too blue to be a normal star, it cannot be 
considered a pulsar counterpart candidate because of the 
large offset (3'.'4) from the radio position, provided the ra- 
dio astrometry is as accurate as claimed by Lommen et 
al. (200C). Like most other blue objects in the field at this 



^ http://www.eso.org/observing/dfo/quality/FORS2/qc/ 
photcoefT/ photcoeffs_fors2.html 



magnitude level (we found at least five objects roughly of 
the same color), it is most likely extragalactic, and simi- 
lar to blue objects found in the HDF images. The object 
marked "6" in the upper right of the Rg image in Fig. ^ 
is only marginally detected at i?s=27.45±0.61. It is not 
detected in the other bands and may well be an artefact 
of the reductions. It is also too far (3'.'1) from the radio 
pulsar position to be an optical counterpart candidate. 

3. Discussion 

Based on the apparent similarity of the radiation proper- 
ties of MSPs and ordinary pulsars, and assuming the same 
mechanisms to generate optical emission in both types of 
pulsar, we expected to detect the optical counterpart of 
PSR J0030+0451 at the m 26 visual magnitude level, 
assuming a simple scaling of E/(P from known optical 
fluxes of ordinary pulsars. Our observations were deeper, 
but still we did not detect any reliable counterpart candi- 
date. 
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Fig. 2. Radio and X-ray observations of PSR J0030+0451 with the A recibo telescope (Lommen et al. ^000[ ), ROSAT 
([Becker et al. 2000|) , and XMM-Newton (|Bccker fc Aschenbach 2002[) , as well as VLT upper limits in the BVRs op- 
tical bands. Solid, dashed, and dot-dashed lines show the best spectral fits of the XMM-Newton data with three 
alternative two component spectral models: BB+PL, broken PL, and BB+BB, respectively (BB stands for black- 
body). Parameters of the fi ts and a mean unabsorbed integrated flux Fx, approximately the same for all the fits 
( Becker &: Aschenbach 2002 ), are indicated in the plot. All fits are acceptable at the same significance level in the 
^ 0.3 — 7 keV range, indicated by a horizontal bar with arrows at the upper-left, and their ~90% uncertainties are 
shown as shaded regions. The unabsorbed model spectra are extrapolated toward the optical range. Previous ROSAT 
data contain no spectral information and provide only an integrated flux indicated in the figure by a bold cross 
( Becker et al. 20do| ), roughly consistent with the more recent XMM-Newton results. The VLT upper limits, marked 
by end bars, show that any nonthermal PL component obtained from the above spectral fits of the X-ray data must 
be strongly suppressed in the optical range. This implies either a cutoff or a strong break in the 0.003 — 0.1 keV range. 
The purely thermal BB+BB model is consistent with our optical upper Hmits. Further details of this plot are discussed 
in Sect. ^. 



To try to understand what our optical non-detection 
implies, we have plotted the available information about 
the multiwavelength spectrum of PSR J0G3G+0451 includ- 
ing radio, optical, and X-ray data in Fig. ^. For the X- 
ray region we have included preliminary results of recent 



XMM-Newton observations ( [Becker fc Aschenbach 2002] ) 
where the pulsar was clearly detected in the 0.3 — 7 keV 
range. The detection range is shown by a horizontal bar 
with arrows in Fig. ^ The overall spectrum compares well 
with multiwavelength spectra of ordinary pulsars (e.g., 
Koptsevich et al. 200l| ), i.e., the pulsar flux is higher in 



the radio and fades toward the X-ray range. This can 
be explained by different emission mechanisms in the ra- 
dio (coherent) and at shorter wavelengths (non-coherent). 
However, a more detailed inspection of the optical/X-ray 
range reveals a feature which has not been seen for or- 
dinary pulsars. For the latter, the optical flux is usually 
close to an extrapolation of the nonthermal high energy 
tail of the X-ray emission, usually described by a power- 
law (PL). From Fig. || it is clear that this is not the case 
for the PSR J0030+0451. 
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Fig. g shows that in the XMM-Newton range the 
data can be fitted equally well by three different two- 
component spectral models: a blackbody + power-law 
model (BB-I-PL), a broken (or curved) PL model, or 
a model based on two different blackbodies (BB-I-BB) 
with TVh < 2.5 X 10^° cm~^. From the X-ray data 
alone it is difficult to discriminate between the three 
models, although the sharp X-ray pulse profile perhaps 
favors the domination of a nonthermal PL component 



( [Becker k Aschenbach 20021 ). The VLT upper limits al- 
low additional constraints using an extrapolation of the 
X-ray model spectra toward the optical range. 

From Fig. |2| it is seen that the low-energy extension 
of the BB-I-PL fit overshoots the optical flux upper lim- 
its of PSR J0030-f0451 by -5 orders of magnitude. The 
broken PL extension, which implies a flatter spectrum in 
the soft X-ray energy band, is also 3-4 orders of magni- 
tude higher. We have not corrected the data for inter- 
stellar extinction, which is low and does not play any 
significant role at such large differences. The strong sup- 
pression of the PL components in the optical range sug- 
gests that these two models either should be ruled out, 
or that their PL components must have a strong break 
or even a cutoff at a photon energy somewhere in the 
0.003 — 0.1 keV range. Such a situation has never been 
observed for any ordinary pulsar. For example, the spec- 
tral index (defined as Fi, cx v~"'") of the young Crab 
pulsar changes from ~ 0.5 for soft X-rays to zero in the 
FUV/optical/ncar-IR range (e.g., ^oUerman ct al . 200C; 



SoUerman & Flyckt 2002). For the relatively young Vela- 



pulsar (Mignani & Caraveo 2001), the middle-aged PSR 
B0656+14 ([Koptsevich ct al. 2001D and PSR B1055-52 
(Pavlov ct al. 2002), and even for the old ordinary pul- 
sars PSR B0950+08 ([Zharikov et al. 2002D and PSR 
B1929+10 (iMignani ct al. 2002aD , the extension of the PL 
X-ray component matches the optical flux, suggesting the 
same mechanism of nonthermal emission in the optical 
and X-ray ranges. It is not clear what could be the reason 
for such a strong change in the nonthermal spectral slope 
of PSR J0030-I-0451 from negative in X-rays to positive 
in the optical range. As seen from Fig. ^, the X-ray PL 
fits and our upper limits exclude any flat extension toward 
the optical range, even if we place a break point energy at 
the lower boundary of the XMM-Newton range ~ 0.1 keV. 
Note that the low-frequency extensions of the PL compo- 
nents overshoot even the radio fluxes. This is also not typ- 
ical for ordinary pulsars. One can assume that nonthermal 
emission is due to synchrotron radiation of relativistic par- 
ticles in the magnetosphere of the pulsar. In this case we 
obtain for the simplest monochromatic particle distribu- 
tion over the energy a spectral flux of F^, cx v^^^ in the low 
frequency range below a maximum frequency t'm oc Bj^. 
Here B is the magnetic field, and 7 is the gamma-factor of 
the emitting particles. From the spectral shape suggested 
by the X-ray data and optical upper limits it is natural to 
put near the maximum of the X-ray spectral flux at 
the low boundary of the XMM-Newton range. At typical 
gamma- factors of primary and secondary relativistic parti- 



cles in pulsar magnetospheres, i.e. ~ 10^ and ^ 10, respec- 
tively, and for the period of PSR J0030-f 0451 - 4.9 ms, 
the same peak frequency value is predicted by a model of 
synchrotron emission from the pulsar light cylinder sug- 
gested by Malov (2001). For the expected synchrotron flux 
below the adopted j/m value (see Fig. ^ dotted lines) we 
would likely hint the optical counterpart, but we did not. 

The purely thermal BB-I-BB spectral model is consis- 
tent with our upper limits without any additional assump- 
tions. Its Rayleigh- Jeans tail is about 6 stellar magnitudes 
fainter in the optical than our upper limits, and would 
hardly be detectable with present telescopes. Thermal 
photons can be emitted by hot polar caps of the pulsar. 
The two-blackbody fit indicates a non-uniformity in the 
temperature distribution over the caps. This could be de- 
scribed by heat propagation over the surface of the neu- 
tron star out of a hot cap core, including neutron star at- 
mosphere effe cts, as has been do ne in the case of the MSP 
J0437-4715 ( ^avhn ct al. 2002[ ). An additional faint PL 
component is required to fit an excess over the thermal 
emission at high energy X-rays from PSR J0437— 4715. 
If the same would be true for PSR J0030-I-0451, it could 
be brighter in the optical than estimated from the simple 
BB-I-BB model due to a contribution from a similar non- 
thermal component of magnetospheric origin. Deeper X- 
ray observations are probably needed to detect this com- 
ponent in the high energy tail of the PSR J0030-I-0451 
spectrum. The similarity of X-ray and radio pulse profiles 
of this pulsar suggests that radio and X-ray peaks are in 



phase (Becker & Aschenbach 2002), although direct tim- 
ing to confirm this has not yet been done. In the frame 
work of the thermal model this means that radio emis- 
sion is generated close to the polar cap surface and the 
similarity of the pulse shapes is likely caused by the same 
geometry of the emitting regions. 

Based on our upper limits we can constrain also the ef- 
ficiency of converting spindown power of PSR J0030-f 0451 
to optical emission. The luminosity in the B band is 
Lb = incfFBAvB ^ 6.3 x 10^^ d^go erg s"\ and hence 
the optical efficiency tjb — Lb/E ^ 1.9 x 10^^c?230- 
Here ^230 — c?/230 pc is the normalized distance to 
PSR J0030-K0451. This upper hmit is about half the effi- 
ciency of the middle-aged PSR B0656-H14 (at 500 pc) and 
exceeds the efficiencies of the Geminga and Vela pulsars 
by about 1 and 2 orders of magnitude, respectively (see, 
e.g., Zharikov ct al. 2002 ). It is interesting to note that 
the expected efficiency in the BB+BB model in Fig. ^ is 
2-3 orders of magnitude lower than the upper limit derived 
above from our optical data, and that the efficiency in the 
BB+BB model is comparable to that of the Vela pulsar. 
The comparison of the efficiencies makes sense here, since 
the thermal emission from hot polar caps and the non- 
thermal magnetospheric emission of the Vela-pulsar, both 
are powered by relativistic particles produced in magne- 
tospheres of rapidly rotating NSs. Thus, we see that the 
optical efficiency of the PSR J0030-I-0451, as derived from 
our and X-ray observations, is not unusually low, but is 
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compatible with the efficiency range of ordinary pulsars 
detected in the optical band. 

It has been shown for ordinary pulsars that 
spectral index appears to became steeper with pul- 
sar age in the optical range (Koptsevich et al. 2001 



Mignani & Caraveo 2001), while it flattens in gamma rays 
(e.g., ^hearer fc Golden 200^ ). It has been noticed also 
across a restricted set of young and middle-aged pul- 
sars detected in the optical and gamma regions that the 
gamma-ray efficiency increases with age, while the r everse 
is true for the optical efficiency ( Goldoni et al. 1995 ) . This 
would suggest that there is a reprocessing of the gamma- 
photons into the optical in pulsar magnetospheres, and 
that it is more efficient for younger pulsars than for older 
ones (e.g., Shearer fc Golden 2002| ). Thus, it would be 
not surprising that very old PSR J0030+0451 is fainter 
in the optical than we expected. However, recent opti- 
cal studies of old ord inary pulsars (^harikov et al. 2002 ; 
Mignani et al. 2002a; Mignani et al. 2002b) have revealed 
nonmonotonous behavior of the optical efficiency vs age 
with a minimum at t 10* — 10^ yr and further increase 
towards higher ages r ^ 10^ yr. Old pulsars can be ac- 
tually much more efficient than the middle-aged ones and 
produce the optical photons with almost the same effi- 
ciency as young and energetic Crab-like pulsars. In this 
context low optical brightness of the MSP J0030-h0451 
remains puzzling. 

A clue to what dominates the X-ray and optical 
spectrum would hopefully be found from observations of 
PSR J0030+0451 and other MSPs in the FUV and, in par- 
ticular, in the EUV range. Even deep upper limits in these 
ranges would help to understand how strongly the mul- 
tiwavelength emission and radiation properties of MSPs 
differ from those of ordinary pulsars. 
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